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Abstract 

The Manickam-Miklos-Singhi Conjecture states that when n > Ak, every multi- 
set of n real numbers with nonnegative total sum has at least (^Zi) A;-subsets with 
T-H nonnegative sum. We develop a branch-and-cut strategy using a linear programming 

formulation to show that verifying the conjecture for fixed values of /c is a finite prob- 
lem. To improve our search, we develop a zero-error randomized propagation algorithm. 
Using implementations of these algorithms, we verify a stronger form of the conjecture 
^ for ah < 7. 

> 



£j 1 Introduction 

^ Given a sequence xi, . . . ,x„ of real numbers with nonnegative sum ^^^=1 natural to 

\Q ask how many partial sums must be nonnegative. Bier and Manickam [4] showed that every 

^ nonnegative sum X^iLi least 2"~^ nonnegative partial sums. Manickam, Miklos, 

^ and Singhi [10, 11] considered the situation when each partial sum has exactly k terms (we 
call such partial sums k-sums), and they conjectured there are many nonnegative /c-sums 

O when n is sufficiently large. 

Conjecture 1 (Manickam-Miklos-Singhi [10, 11]). Let k > 2 andn > Ak. For all xi, . . . ,x„ G 
M such that Yl^=i^i — 0; there are at least subsets S C [n] with \S\ = k such that 

The bound {^ZD is sharp since the example xi = n — 1 and Xi = —1 for i I has sum 
zero and a k-sum is nonnegative exactly when it contains xi. The bound n > 4A; is not 
necessarily sharp. 

In this paper, we develop a computational method to solve fixed cases of the conjecture 
based on a poset of /c-sums and a linear programming instance. In particular, we prove a 
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stronger statement than the conjecture for all k < 7. This leads us to formulate a stronger 
form of the conjecture (see Conjecture 2). 

Conjecture 1 is trivial for k = 1 and is a simple exercise for k = 2. Marino and Chi- 
aselotti [12] proved the conjecture for k = 3. Bier and Manickam [4] showed there exists 
a minimum integer f{k) such that for all n > f{k), there are at least (^I^) nonnegative 
fc-sums in any nonnegative sum of n terms. Subsequent results [3, 4, 5, 6, 15] give several 
exponential upper bounds on f{k). Alon, Huang, and Sudakov [1] showed the polynomial 
bound f{k) < mm{33k^ , 2k^} . Chowdhury [7] proved /(3) = 11 and /(4) < 24. 

We prove /(4) = 14, /(5) = 17, /(6) = 20, and /(7) = 23. For A; < 7 and n < f{k), 
we find that the nonnegative sums with the fewest nonnegative /c-sums have the following 
structure: for positive integers a and b summing to n, let Xi = ■ ■ ■ = x;, = a and x^+i = 
■ ■ ■ = Xn = —b. When 3k < n < f{k), the extremal examples have a = 3, 6 = n — 3, and 
("'^^) nonnegative fc-sums. This leads us to make the following conjecture. 

Conjecture 2. Let be the smallest integer such that > (^'li^) ■ Then f{k) = N^, 

and hence^ 

lim = lim ^ ^ 3.147899. 

In the next section we develop some necessary notation and preliminary results. In Sec- 
tion 2, we place a symmetry-breaking constraint on the vectors and investigate a natural 
poset among the fc-sums. In Section 3, we develop our branch-and-cut strategy for verify- 
ing the conjecture. We then modify this algorithm by adding a randomized propagation 
algorithm in Section 4. In Section 5 we discuss stronger forms of the conjecture. 



1.1 Preliminaries 

Let ai, . . . , be real numbers and ei, . . . , be positive integers. We denote by al^ ■ ■ ■ 'h^ 
the vector in with the first Ci coordinates equal to ai, the next 62 coordinates equal to 
02, and so on until the last coordinates are equal to a^. Using this notation, the vector 
{n — 1)^ (—1)""^ is conjectured to achieve the minimum number of nonnegative /c-sets when 
n > 4k. 

For a vector x = (xi, . . . , x„) G and a k-set S, define 0"5(x) = Y2ies "^^^ number 

of nonnegative fc-sums of x is Sfc(x) = G (j-^) : o"5(x) > o| . For small values of k and n, 

we will compute the minimum integer g{n, k) such that all vectors x G M" with nonnegative 
sum have Sfc(x) > g{n,k). By the sharpness example, g{n,k) < always. Define the 

deficiency function as g{n, k) = (Izl) " di^i 

Since removing the least coordinate from a vector x G M" results in an (n — l)-coordinate 
vector with nonnegative sum, the function 5f(n, k) is nondecreasing in n. Bier and Man- 
ickam [4] proved that g{n,k) = when k divides n by using a theorem of Baranyai [2]. 



^This limit was computed by first simplifying — {^_\) = to (a; — k){x — k — l){x — k — 2) — k{x - 

l){x — 2) =0 and taking the real root for x of the resulting cubic. 
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The following lemma of Cliowdhury [7] reduces the problem of computing f{k) to verifying 
g{n, k) = for a finite number of values n. 

Lemma 3 (Chowdhury [7]). If g{n, k) = 0, then g{n + k,k) = 0. 

This lemma allows the computation of f{k) by first computing a finite number of values 
of g{n, k). We develop a finite process to verify g{n, k) > t for fixed values of k and n, and 
hence by Lemma 3 we determine f{k) when we find that g{n,k) > (JtI) for k consecutive 
values of n. Our computation to verify g{n, k) > t for fixed integers n, k, and t is designed 
as a search for a vector x e with Sfc(x) < t. By constraining the order of the coordinates 
in a vector x G M", we are able to make significant deductions about the sign of the fc-sums 
of X, which makes the computation very efficient. We first discuss this constraint on the 
coordinates of x in Section 2, and then discuss the algorithm in Section 3. 



2 The Shift Poset 

To better control the vectors x e with nonnegative sum, we make a simple ordering con- 
straint that adds significant structure. Let Fn be the set of vectors x G M" with nonnegative 
sum and with the sequence Xi nondecreasing. That is, 

n 




A — 1 



i=l 

Denote by [n] the set {1, . . . ,n} and by ('^') the collection of fc-subsets of [n]. For two 
k-sets S = {ii < 12 < ■ ■ ■ < ik} and T = {ji < 32 < ■ ■ ■ < jfc}, let 5" ^ T if and only if < ji 
for all £ G {!,..., k}. We call this poset over ('^^) the shift poset. Since we write Xi, . . . ,Xn 
from left to right in nondecreasing order and associate a k-set S as with the values Xi for 
i G S*, then when S* ^ T we say S is to the left of T and T is to the right of 5*. Observe that 
a relation ^ T is a cover relation if and only if 5* \ T = {i}, T \ S = {j}, and i = j — 1. 

Since we restrict vectors in F„ to have nondecreasing values as the index increases, S 
imphes o'si'x.) > cttIx) for all x G F„. Therefore, if crT(x) is nonnegative, then so is asi'x). 
Similarly, if o"5(x) is strictly negative, then so is crT(x). 

Given S G {^f), let the left shift of S be CkiS) = {T E {^l^) : T y S} and the right shift 
of S be TZkiS) = {T G ('J) : S^T}. Given a set A C (W), the left closure of A, denoted 
Ck{A), is the union of all families Ck{S) over all sets S E A. Similarly, the right closure of 
A, denoted 7lk{A), is the union of all families 7lk{S) over all sets S E A. 

The shift poset has many interesting properties, including the fact that it is a lattice, 
which will be critical in later calculations. Fix any list J-" = {5*1, . . . , Si} of fc-sets where 
Sj = {ij^i < ■ ■ ■ < ij^k}- The meet of J-', denoted AJ-" or A^^^S^, is the fc-set T such that for 
all T' where Sj >z T' for all j G [i], then T ^ T'. The join of J-", denoted VJ-" or A^.^Sj, is 



the fc-set T such that for all T' where T' >z Sj for all j G [i], then T' >z T. The meet and 
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join can be computed as 

A-F = A'j^.Sj = { max{i,-t : j E [£]} : t e [k]}, 
= yU^, = { min{z,-i :je[i]}:te[k]}. 

Observe that the fc-sets to the left of all 5*1, ... , Se are exactly the fc-sets to the left of the 
join V^^j^S'j, and the fc-sets to the right of all S*!, ... , 5"^ are exactly the k-sets to the right of 
the meet A^^^Sj. That is, 

e I 

i=i i=i 

The shift poset preserves order with two standard linear orders of fc-sets: the lexicographic 
and colexicographic orders. Let lex(S') and colex(S') be the lexicographic and colexicographic 
rank, respectively, of a fc-set S* = {zi < ■ ■ ■ < 4} ^ {1, • • • ,n}, defined as 

where io = by convention. These functions are bijections from ('^') to the set {0, . . . , (^) — 1} 
with the property that when S we also have lex(S') < lex(T) and colex(5') < colex(T). 
In addition to these ranking functions, the lexicographic and colexicographic orders also 
have efficient successor and predecessor calculations. We use these methods to iterate over 
all fc-sets to find ^-maximal or ^-minimal elements. 

Finally, we note that the shift poset is isomorphic to the dominance poset over composi- 
tions of n + 1 into k + 1 positive parts. Given two compositions of n + 1 into k + 1 positive 
parts, where A = {ai, . . . , a^+i} and B = {bi, . . . , bk+i}, A < B in the dominance order 
when ^l^itti < Yli=i^i fo'^ t G [k]. Note that the dominance order, also called the 
majorization order, is usually defined over partitions, but the defining inequalities are the 
same for compositions (see [14]). 

For S = {ii < ■ ■ ■ < ik}, let As = {ai, . . . , Ofc+i} where ai = ii, ag = i^ — ie-i for 
2 < £ < k, and a^j^i = n + 1 — ik. The image As has the property that Yll^i a£ = it for 
t < k and Yl^=i o,£ = n + 1. Therefore, if 5" ^ T in the shift poset, then As < At in the 
dominance order. 

2.1 Counting Shifts 

The functions Lk{S) = \Ck{S)\ and Rk{S) = |7^fc(S')| count the size of the left and right 
shifts, respectively. These shift functions are quickly computable using a recursive formula 
due to Chowdhury [7] which we prove for completeness. If S" = {ii, . . . ,ik} and j < k, let 
Sj = {ii, . . . , ij}. We define Lj{S) = \Cj{Sj)\. 
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Proposition 4 (Chowdhury [7]). For a k-set S = {ii < i2 < ■ ■ ■ < ik} ^ {^f), 

^'(-)K':)-C'rO-|:^'<^'(V-T)- 

Proof. Observe that £^(5) C (l^^^l). We will count the sets T = {ji < ■ ■ • < in {^f) \ 

Ck{S) as ('*^*^) + Yl^=i ^ii^^){^''k-e^) ■ Every such set T has some minimum parameter 
i > such that j^+i > ie+i- If i = 0, then the minimum of T is strictly larger than ii and 
there are exactly (*'°^*^) such sets. When i > 1, we have jt < «t for all t < £. These sets 
can be constructed from the Li{Se) ^-subsets of [n] that are to the left of Se and extended 
by exactly (^^^l^') ~ £)-subsets of {ie+i + 1, . . . , 4}- ■ 

Observe that if S* = {zi < ■ ■ ■ < ik}, then by symmetry for T = {n + 1 — ig : i E [A;]} we 
have Rk{S) = L^iT). While computing L^i^S) for all fc-sets, we first compute and store the 
values of Lj{Sj) for 1 < j < /c. 

2.2 Required Negative Sets 

We use the shift function -^^(5) to determine that certain /c-sums must be negative in order 
to avoid t nonnegative fc-sums. 

Observation 5. Let n,k,t be integers. If S E (^^'^) has Lk{S) > t, then every vector yi G 
with Sfe(x) < t has (x) < 0. 

Lemma 6. Fix integers n,k,t with t < (^'I^J) and a k-set S G ('^') with 1 E S. If Lk{S) + 
g{n — k,k) > t, then every vector x G F„ with Sfc(x) < t has ^^(x) < 0. 

Proof. Fix a vector x G -F„ with Sfc(x) < t. The fc-set T = {l,n — k + 2, . . . ,n} has 
Lk(T) = Qlj), so by Observation 5, o"t(x) < 0. Since Y17=i^i — 0' have '^^=2'^^ ~ 
Tlii^^T^i — 0- ^ family of nonnegative fc-sums over X2, ■ ■ ■ , Xn-k- By the definition 

of g, \A\ > g{n — k,k). For all Q E A, the minimum element of Q is at least 2 and hence Q 
is not to the left of S, since the minimum element of 5* is 1. Therefore, if crsi'^) ^ 0, then 
■Sfc(x) > \Ck{S) \ + \ A\ > Lk{S) + g{n — k, k) > t, a contradiction. ■ 

Similar inferences, discussed in Section 3.3, will be important to the performance of our 
algorithms. 

3 The Branch-and-Cut Method 

Our method to verify the conjecture centers on searching for vectors x G which are 
counterexamples (i.e. Sfc(x) < (^I^))- In addition, for values of n less than /(/c), we will 
search for vectors with the fewest number of nonnegative fc-sets. For integers n, k, and t, 
we say a vector x G is (n, k,t)-bad if Sfc(x) < t. Using the shift poset, we will determine 
properties of such an (n, k, t)-bad vector and use that to guide our search. 
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For a vector x e F„, the /c-subsets of [n] are partitioned into two parts by whether the 

associated fc-sums arc nonnegative or strictly negative. Let contain the k-sets S with 
o's(x) > and contain the the k-sets T with (Tt(x) < 0. Observe = Ck{C^) and 
= 7lk{C~). Our computational method focuses on testing whether vectors x G F„ exist, 
given certain contraints on C+ and C~. 

Definition 7. Given families and A~ of k-sets, let the linear program P{k,n, A'^ , A~) 
be defined as 



Vi e {!,..., n-1} 
ys e A+ 

WeA- 



A vector x e F„ has D and C~ D TZk{A~) if and only if an appropriate scalar 

multiple of X is a feasible solution to P{n, k, A'^, A~). 

Due to the following observations, we can verify that .Sfc(x) > t for the infinite set of 
vectors x in F„ by searching for partitions of ('^') that correspond to the nonnegative and 
negative fc-sums of x. 

Observations. //£fc(^+)U7^fc(^-) = ('J), then every feasible solution x to P{n, k, A'^, A ) 
has sa;(x) = In particular, a feasible solution to P{n, k, A'^ , A~) is {n, k,t)-bad if 

and only if \Ck(yA'^)\ < t. 

Observation 9. If P{n, k, A'^ , A~) is infeasible, then no vectors e F„ has D Ck{A'^) 
and C~ 3 7lf,[A^). 

Observation 10. If\Ck{A'^)\ > t, then Sfe(x) > t for all feasible solutions to P{n, k, A'^ , A~). 

By searching for sets A^,A~ such that a solution to P{n,k,A'^,A~) has Sjk(x) < t, we 
either find these vectors or determine that none exist. Thus, determining g{n, k) is a finite 
problem, and by Lemma 3 determining f{k) is a finite problem. 

3.1 The Search Strategy 

Fix n, k, and t < i^Z}^- We will search for sets A'^ and A~ such that the solution x to 

P{n,k,A-^,A-) has Sfc(x) < t. If U7^fe(^-) is a partition of (^J) with |A(^^)| < 

t, then every feasible solution to P{n,k,A'^iA~) is (n, /c, i)-bad. The reason we use the 



P{k,n,A\A ): minimize x\ 

n 

subject to '^^Xi > 

i=l 

ies 

Y^xi < -1 

xi, . . . ,Xn eR 
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objective of minimizing xi in the linear program is to try and distance the optimal solution 

from the conjectured sharp example of Xi being a large positive value and X2, ■ ■ ■ , x„ being 
small negative values. In practice, when an (n, k, t)-bad vector exists we discover it before 
Ck{A^) and Tlk^A-) partition ('^1). 

During the algorithm, we will store two collections of /c-sets, and B~ , called branch 
sets. We initiahze B'^ — B" — and always B^ C and B~ C A". The difference 
between B~ and A~ , for instance, is that the sets in B~ are chosen to have negative sum, 
but the sets in A~ \ B~ are sets that must have negative sum for all (n, k, t)-bad vectors x 
with C+ D i3+ and C" 3 B". 

The search procedure BranchAndCut (Algorithm 1) is recursive, taking parameters 
n,k,t,B+,B-,C* under the requirement that >Cfe(B+) U TZkiB') U C* is a partition of ('J). 
The collection C* = ('^1) \ (£fc(-B+) U TZkiB-)) contains the A;-sets S where (T5(x) is not 
immediately decided by the constraints on B^ and B~ . 

We will refer to each recursive call to BranchAndCut as a search node, and at each 
search node we perform the following three actions: 

1. Determine the ^-maximal sets S' G C* such that all vectors x with C+ 2 B"*" U {S} 
have sa;(x) > t and add these sets to A~ . 

2. Test that the linear program P{n, k, B"*", A~) is feasible (prune if infeasible). 

3. Select a set 5" e C* and branch on the choice of placing S in B^ or B~ , creating two 
new search nodes. 

One crucial step for the first action performed at each node is computing the number of 
/c-sets in C* that are also to the left of a set 5" or to the right of a set S. Define the functions 
Ll{S) = \Ck{S) \ Ck{A^)\ and Rl{S) = [R^iS) \ 7^,(^-)|. 

The algorithm SelectBranchSet is a method for selecting a set G C* for the branch- 
ing step is called a branching rule. A good branching rule can significantly reduce the number 
of visited search nodes. However, it is difficult to predict which set is the best choice. Any 
branching rule that selects a A;-set S from C* will provide a correct algorithm, but these 
choices can greatly change the size and structure of the search tree. Based on our experi- 
ments and choices of several branching rules, we found selecting a /c-set S that maximizes 
min{L^(S'), -R^(S')} was most effective. This rule ensured that both branches removed as 
many sets from C* as possible. 

In Section 3.3, we define the algorithm PropagateNegative (Algorithm 2). Before 
that, we must discuss how to efficiently compute L*^.{S) and Rl{S). 

3.2 Computing Intersections with C* 

For this discussion, fix two famihes ^+ and A~ with C* = (^J) \ (£fc(v4+) MkiA')). We will 
use three methods to compute \Ck{S) r\C*\ and \TZk{S) nC*|: breadth-first-search, inclusion- 
exclusion, and iterative updates. In all techniques, we will store markers for whether a set 
T is in Ck{A+), TZkiA-), or C*. 
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Algorithm 1 The recursive algorithm Branch ANDCl]T{n, k, t, B~^, ,C*)- 
Require: jCkiB+), TZkiB-), and C* partition (IJ). 

Ensure: Returns an (n, k, i)-bad vector x with C+ D B^ and C~ D B~ if it exists, 
if \C{B+)\ > t then 

return Null 
end if 

A'^,A'',C* ^ PROPAGATENEGATiVE(n, k, t, B+, B~ ,C*) {See Algorithm 2) 
Find optimal solution x to k, i A~) 
if X = Null or Sfc(x) < t then 

The linear program is infeasible or found a counterexample. 

return x 
end if 

S ^ SelectBranchSet(^+,^-,C*) 
Bt^B+, B^^B-U{S}, Cl^C*\nk{S) 
X BRANCHANDCUT(n, /c, t, iSj*", Bl.Cl) 
if X = Null then 
Bt^B+U{S}, B^^B-, C*,^C*\Ck{S) 

X BRANCHANDCUT(n, k, t, Bj, B2 , C2) 

end if 
return x 



Breadth-First-Search. When k is small, it is reasonable to store the Hasse digraph (an 
edge 5" ■<— T exists if 5* ^ T is a cover relation) of the shift poset in memory and use 
breadth-first-search to count the number of sets to the left of -S" which are in C*. 

Inclusion-Exclusion. We can use the lattice structure of the shift poset to compute the 
functions Ll{S) and Rl{S). First, we compute the values of Ll{S) with increasing colex 
rank and values of Rl{S) with decreasing colcx rank. With this order, we have access to 
Ll{T) whenever T ^ S or Rl{T) whenever S ^ T. Let J^iiS) be the collection of sets T 
that cover S, and let J^r{S) be the collection of sets T that S covers. Then, 

$^ ACTUS) 

When k < 4, the breadth-first-search strategy is faster than the inclusion-exclusion strat- 
egy. However, for A: > 6, the inclusion-exclusion strategy is the only method that is tractable. 
For k = 5, the two strategies are too similar to demonstrate a clear preference. 

As k grows, the number of terms of the inclusion-exclusion sum grows exponentially, 
leading to a large amount of computation required for every set S & C*. Our next method 
computes the values of Ll(S) when a small change has been made to B'^ using a smaller 
amount of computation for every set S & C*. 
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Updates with B'^. Observe that sets are added to only when placing a set into B+. 

Therefore, there arc many fewer sets in A'^ than in A~ . Further, if the most recent branch 
selected a A;-sum to be negative, then the values L^S) did not change for any sets still in 
C*. Therefore, we need only update the values of Ll{T) for T E C* after branching on a set 
S and placing S in For all such T, observe that 

LKT) ^ Ll{T) - Ll{T y S) 

properly updates the value of Ll{T) when adding the set S to B'^. Specifically, the sets that 
are removed from Ck(T) \ ^{A'^) to Ck{T) \ Ck{A'^ U {S*}) are exactly those to the left of 
both T and S but not to the left of any elements in A'^ . Such /c-sets are exactly those in 
Ck{T y S)\ which are counted by L*(T V S). 

This method is not efficient for computing Rl{S), since too many sets are being added to 
A~ during the propagation step (see the next section). Thus, we use this method to update 
Ll{S) only when a single set has been added to B'^. If we need to recompute all values of 
Ll.{S) after a significant change to B~^, we use the inclusion-exclusion method. 

3.3 Propagation 

Given a set of branch sets B^ and B~ , we want to determine which sets S E C* must have 
(T5(x) < for all (n, A;,t)-bad vectors x which are feasible in P{n, k, A'^ , A~). 

Recall that by Observation 5 and Lemma 6, any set S with Lk{S) > t, or 1 E S and 
Lk{S) + g{n — k,k) > t, has (^^(x) < for any (n, k, t)-bad vector x. This applies regardless 
of the previous choices of sets placed in B'^. When is non-empty and we have access to 
the function Ll{S), we may be able to find more sets where asi'x) < for any (n, k, i)-bad 
vector x with D B^. 

Lemma 11. If S is a k-set with Ll{S) + \jCk{B'^)\ > t, then all feasible solutions x to 
P{n, k, B+ U {S}, B-) have Sjk(x) > t. 

Proof. Observe that such vectors x have C+ D U Ck{S) and |>Cfc(6+) U Ck{S)\ = 

+ \Ck{S) \ = \Ck{B+)\ + Ll{S) > t. 

The PropagateNegative algorithm (Algorithm 2) iterates on all sets S in C* using 
lexicographic order and whenever a set S satisfies the hypotheses of Lemma 6 or Lemma 11 

the set S is added to A~ . Since all sets T ^ S have T <icx 5", the set T was considered 
before S and did not satisfy the condition for addition to A~ . Hence, we construct A~ as a 
^-maximal set (other than possible comparisons within B~). 

3.4 Results Using BranchAndCut 

Using BranchAndCut (Algorithm 1), we computed g{n, k) for k e {4, 5, 6} and k < n < 
4/c + 1, and we verified that f{k) = 3k + 2 ioi k e {4, 5, 6}. 

Theorem 12. /(4) = 14, /(5) - 17, and /(6) = 20. 
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Algorithm 2 PROPAGATENEGATiVE(n, /c, i, B+, C*) 
Require: 7^fc(B-), and C* partition (1^1) and \jCkiB+)\ < t. 

Ensure: Returns B~^, Aq , and Cq such that any (n,A;,t)-bad vector x with D B'^ and 
C- D B- also has C-DA^, and Q = ('^1) \ {Ck{B+) U UkiA^)). 

Ao^B-, Q^C* 

for all sets S E Cq (in lex order) do 

if 1 G 5 and Lfe(S) + g{n -k,k)>t then 
Ao ^ Ao U {S} 

end if 

if Ll{S) + \Ck{B+)\>t then 
A" ^ A U {^} 

end if 
end for 

return . A7. . C,', 



Proof outline. By Lemma 3, we must only verify the values n where f{k) <n< f{k) + k — 1. 
For all n < 3A; + 1, we tested all sums a + h = n and found Sfc(x) for the vector xi = ■ ■ ■ = 
Xb = a, Xh+i = ■ ■ ■ = Xn = —b. Wc let t be the minimum such value Sfc(x) and executed 
BRANCHANDCuT(n, k,t, 0, 0, ('^')), which found no vector with fewer than t nonnegative 
/c-sums. By executing these algorithms in increasing value of n, we have access to g{n — k, k) 
during execution of PROPAGATENEGATlVE(n, /c, B+, B~,C*) (Algorithm 2) when testing 
g{n, k) >tfor n > 2k. Finally, forn G {3A; + 2, . . . , 4A; + 1}, we find g{n, k) > glj) by 

executing BRANCHANDCuT(n, k, Q^J), 0, 0, (^^1)). ■ 

For k — 7, BranchAndCut succeeded in computing g{n, k) for n < 23, but only after 
resorting to parallel computation. Previous computations showed an order of magnitude 
jump in computation time between n = f{k) — 1 and n — f{k), so using BRANCH AndCut 
to verity g{2A, 7) = seemed intractable. 

In the next section, we develop a zero-error randomized propagation algorithm which 
improves the performance enough for us to compute /(7). 

4 Propagation With Randomness 

When branching in our branch-and-cut method, we select a set S in C* and make a tem- 
porary decision of o's(x) < — 1 or o's(x) > 0. Then during our propagation step, the 
PropagateNegative algorithm (Algorithm 2) places as many left-most sets into A~ as 
possible using Lemmas 6 and 11. We now describe a condition that allows us to add sets 
into A'^ without branching. 
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Algorithm 3 PROPAGATEPosiTlVE(n, A;, C*) 

Require: 7^fe(^-), and C* partition (^J) and \Ck{A'^)\ < t. 

Ensure: Returns Aq, Aq , and Cg such that any {n,k,t)-had vector x with D A^ and 
C- D A- also has C+ D A^, A^ =A-, and C* = ('J) \ {Ck{A^) U UkiA^)). 

Aq,Ao,Co ^ PROPAGATENEGATWE{n,k,t,A'^,A~,C*) 
updated <— True 

while updated do 

update ■<— False 

for all sets S E Cq (in reverse lex order) do 
if V{n, k, A+,A-U {S}) is infeasible then 
update ^ True 
A^ ^ At u {S}. 
C*^C*\C,{S) 

Ao,Ao,Co PROPAGATENEGATIVE(n, k, t, Aq , Aq , Co) 
end if 
end for 
end while 
return Aq,Aq,Cq 



Given a set 5" in C*, we can test the hnear program for feasibility when S is added to A 

or 

Observation 13. If S is a k-set where P{n, k, A'^ , A^ U {S}) is infeasible, then all vec- 
tors that are feasible solutions to P{n, k, A'^, A^) are also feasible solutions to P{n, k, A'^ U 
{S},A-). 

Observation 14. If S is a k-set where P{n, k,A'^ U {S},A~) is infeasible, then all vectors 
that are feasible solutions to P{n, k, A^jA") are also feasible solutions to P{n, k, A'^, A~ U 
{S}). 

Given these two lemmas, we could test every set in and determine which sets should 
be added to A'^ or A~ . The PropagatePositive algorithm (Algorithm 3) places as many 
right-most sets into A'^ as possible using Lemma 13. 

Since PropagatePositive optimizes a linear program for every set in C*, this algorithm 
is too slow to be effective within BranchAndCut. However, since PropagatePositive 
includes a call to PropagateNegative after every addition to B+, it is possible that a 
single call to PropagatePositive results with the linear program V{n,k,AQ,AQ) infea- 
sible. Since PropagateNegative adds negative sets to A^ ^ there are likely more sets 
5* e C* where adding S to A^ leads to an infeasible linear program. Such sets are added to 
A'^ by PropagatePositive, leading to more sets which satisfy the conditions in Propa- 
gateNegative. 

Using PropagatePositive in place of PropagateNegative in the propagation step 
of BranchAndCut is slow for even k — b. However, ioi k — 2> and k — 4, using Propa- 
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Algorithm 4 STOCHASTlcPROPAGATlON(n, A;, i, C*) — Randomly test sets for in- 

clusion in A~ and A'^. 

Require: A; > 3, n > k, t < {^zl), and B+,B- C {^l}). 

Ensure: Sets added to A'^ or A^ satisfy the hypotheses of Observations 13 or 14. 
loop 

A~,C* PROPAGATENEGATiVE(n, k, t, A'X*) {see Algorithm 2) 
updated False 
while not updated do 

Randomly select a set S* G C* 
if V{n, k, A' U {S}) is infeasible then 
A+ ^A+U {S}, C* ^C*\ Ck{S), updated ^ True 
if > t then 

return A'^,A~,C* 
end if 

else if P(n, k, A'^ U {S}, A~) is infeasible then 

A-^A-U{S}, C* ^C*\nk{S), updated^ True 
end if 
end while 
end loop 

return A'^,A~,C* 



GATePositive terminates in at most three iterations of the loop and requires no branching. 
The log of this computation is small enough that we present the computation as proofs that 
/(3) < 11 and /(4) < 14 in Appendices B and C. 

While PropagatePositive is unreasonably slow, PropagateNegative is very fast 
to test since we have very quick methods for computing Ll{S). If we can find just one set 
S where adding S to A~ creates an infeasible linear program, then we can add S to and 
likely the next iteration of PropagateNegative will add more negative sets. 

Instead of carefully selecting a set S to peform this test, we simply select a set S in C* 
uniformly at random. By randomly selecting sets and testing the linear program, we may 
quickly find a set that we can guarantee to be in for any (n, k, t)-bad vector x G This 
is the idea for the algorithm StochasticPropagation (Algorithm 4). Simply, we select 
a set from C* at random and test if the set fits the hypotheses of Observations 13 or 14. 
We continue samphng until either (a) we find such a set and add it to A'^ or A", (b) we 
sample a specified number of sets which all fail these conditions, or (c) we reach a specified 
time limit. The sampling limits of (b) and (c) are not listed in Algorithm 4, but are both 
parameters which can be modified in our implementation. 

StochasticPropagation is a zero-error randomized algorithm: it adds sets to A'^ or 
A~ only when previous evidence guarantees that is the correct choice. The only effect of the 
randomness is how many sets actually are determined to be placed in A'^ or A~. 

What is particularly important is that the propagation in STOCHASTICPROPAGATION 
is stronger than PROPAGATENEGATIVE, but it is still slower. For the best performance. 
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we need a balance between the branching procedure and the propagation procedure. This 
balance is found by adjusting the number of random samples between successful updates in 
StochasticPropagation, as well as the total time to spend in that algorithm. Further, 
we can disable the call to StochasticPropagation until a certain number of fc-sets have 
been added to or B~ , thereby strongly constraining the linear program and leading to a 
more effective random sampling process. 

Solving the case k = 7 and n = 22 required more than 100 days of computation time for 
the deterministic branch-and-cut method, but required only a few hours using the random- 
ized algorithm. 

4.1 Results Using StochasticPropagation 

Using BranchAndCut with StochasticPropagation, we computed g{n, 7) for 8 < n < 
29 and verified that f{k) = 3fc + 2 for A; = 7. 

Theorem 15. /(7) = 23. 

The computation for this theorem is the same as in Thoerem 12, except we replace 
PropagateNegative with StochasticPropagation. 

4.2 Implementation 

The BranchAndCut algorithm was implemented in the MMSConjecture project within 
the SearchLib collection^. The software is available for download, reuse, and modification 
under the GNU Public License 3.0. 

For k G {4,5}, we were able to verify the statements using the exact arithmetic solver 
supplied with GLPK [9], but for /c > 6 this method was too slow and instead was verified 
using the floating-point linear programming software CPLEX [8]. The number of search 
nodes for each search, as well as the amount of computation time for each case are presented 
in Tables 1, 2, and 3. Empty cells refer to computations that were too long to complete (but 
only for cases where k divides n) and cells containing " — " refer to experiments reporting less 
than 0.01 seconds of computation time. In addition to the data presented here, all collected 
statistics and sharp examples are available online^. 

Execution times under one day were executed in a single process on a 2.3 GHz Intel 
Core 15 processor. Longer execution times are from parallel execution on the Open Science 
Grid [13] using the University of Nebraska Campus Grid [16]. The nodes available on the 
University of Nebraska Campus Grid consist of Xeon and Opteron processors with a range 
of speed between 2.0 and 2.8 GHz. 

^SearchLib is available at http://www.math.illinois.edu/~stolee/SearchLib/. 
■^All data is available at http://www.math.illinois.edu/~stolee/data.htm. 
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Figure 1: Values of N^/k for 5 < A; < 250 and the line y = limfc_>.oo N^/k. 

5 Sharp Examples, Uniqueness, and the Strengthened 
Conjecture 

In Tables 1, 2, and 3, the right-most column contains descriptions of the vectors x with Sfc(x) 
of minimum value. For n > f{k), the vectors listed had Sfc(x) of minimum value while also 
having crT'(x) < for T = {l,n — k + 2, . . . , n}. Observe that for all n < f{k), the extremal 
examples use only two numbers, and have the form (—6)" where a + b = n. 

Recall Conjecture 2, where we claim f{k) is exactly equal to A^^, where is the minimum 
integer such that (''^'^T^) > {^k~i)- make this conjecture for two reasons. First, the 
example of 3"^^ (— — 3))^ was known to have fewer than {^Zi) nonnegative fc-sums when 
3k < n < Nk by previous authors (for example, see [7] where these vectors were used to 
show a lower bound f{k) > ^), but no examples were previously discovered for n > Nk- 
Second, our search for extremal vectors found no violation to this bound, but also showed 
that the extremal examples have the form (—6)" when n < and k does not divide 
n. We therefore tested all vectors of the form (— for all k < 250. For all numbers n 
with 3k < n < Nk the example 3"~^ (~(^ ^ 3))'^ was the best vector of this form, and for 
Nk < n < Ak the best vector was {n — ly (—1)"^-'^. 

Figure 1 contains a plot of Nk/k and the line y = limk^ooNk/k, to demonstrate the 
conjectured values of f{k)/k. 

Another strengthening of the conjecture that follows from our method (for k < 7) and 
that of Chowdhury (for A; = 3) is that when n > f{k), a vector x G -F„ has Sfc(x) = 
only if xi + Xn-k+2 + ■ ■ ■ + Xn is nonnegative. Essentially, any sharpness example contains 
the same set of nonnegative /c-sums as the sharpness example {n — 1)^ (—1)""^. 

We test this example is essentially unique by searching for an vector x G -F„ where 
(Tt(x) < for T = {1, n — /c + 2, . . . , n} and Sfc(x) < {^Zl)- If ^lo such vector is found, 
then the sharpness example is "unique," since all vectors x G F„ with sa;(x) < have 
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= Ck{T). Define gs{n, k) to be the minimum Sfc(x) over all x G -F„ where cr2-(x) < 0. We 
can use our target value t to find an [n, fc,t)-bad vector x G -F„ with crT(x) < 0. While the 
bound t < {^Zi) ill the hypothesis of Lemma 6 does not hold, observe that since crr(x) < 
the conclusion of the lemma does hold in this scenario. Therefore, PropagateNegative 
(Algorithm 2) remains correct when verifying gs{n,k) > t using a call to BranchAnd- 
CUT(r^,fc,^,0,{T},(M)\7^,(T)). 

Values of gs{n,k) were computed for 3 < < 6, f{k) < n < f{k) + k and are given 
in Table 4. The sharp examples for vectors x G -F„ with o"t(x) < are given in the final 
columns of Tables 1 and 2. Observe that the sharp examples make a "phase transition" 
at n = 4A;, where the sharp examples have the form (—6)° for all n < 4k, but then the 
examples with n > Ak have at least three distinct values. This may be hinting to a deeper 
truth concerning the originally conjectured bound of f{k) < Ak. 

5.1 Conclusions 

We developed two methods for verifying g{n, k) > t, which we used to prove our strength- 
ening of the Manickam-Miklos-Singhi conjecture for all k < 7, extending the previously 
best result of < 3 [7]. Our branch-and-cut method Branch AndCut uses a branching 
procedure which we prove will, in finite time, verify g{n, k) > t or find a vector x G -F„ 
with Sfc(x) < t. Using the randomized propagation algorithm StochasticPropagation, 
we computed f{k) and the values of g{n,k) for all A; < 7 and k < n < f{k) + k. Our 
implementations were not successful in extending our results to larger values of k due to a 
combination of large computation time and memory requirements. 
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A Computation Data 



k 


n 


g{n,k) 


g{n, k) 


Nodes 


GLPK 


CPLEX 


Strong Example 


3 


4 


1 


2 


2 


— 


— 


1' (-3)^ 


3 


5 


3 


3 


2 


— 


— 


32 (-2)3 


3 


6 


10 





8 


— 


— 




3 


7 


10 


5 


2 


— 


— 


25 (-5)2 


3 


8 


16 


5 


2 


— 


— 


53 (-3)^ 


3 


9 


28 





2 


— 


— 




3 


10 


35 


1 


4 


— 


— 


3'^ {-If 


3 


11 


45 





6 


— 


— 


7^ {-4y 


3 


12 


55 





2 


— 


— 


7"" 1-5}' 


3 


13 


66 





2 


— 


— 


4" (-9)* 


4 


5 


1 


3 


2 


— 


— 


1^ (-4)1 


4 


6 


5 


5 


2 


— 


— 


15 (-5)1 


4 


7 


10 


10 


2 


— 


— 


52 (-2)5 


4 


8 


35 





268 


0.47 s 


0.31 s 




4 


9 


35 


21 


4 


— 


— 


2' (-7)2 


4 


10 


70 


14 


28 


0.11 s 


0.03 s 


2« (-8)2 


4 


11 


92 


28 


10 


0.07 s 


0.01 s 


8^ (-3)« 


4 


12 


165 





50 


0.35 s 


0.11 s 




4 


13 


210 


10 


52 


0.82 s 


0.15 s 


31" (-10)3 


4 


14 


286 





30 


0.65 s 


0.10 s 


104 (-4)1" 


4 


15 


364 





24 


0.51 s 


0.10 s 


114 |_4jii 


4 


16 


455 





6 


0.15 s 


0.04 s 


351 31^^ (-16)5 


4 


17 


560 





8 


0.31 s 


0.07 s 


381 4iu (_i3)e 


5 


6 


1 


4 


2 


— 


— 


15 (-5)1 


5 


7 


6 


9 


2 


— 


— 


1« (-6)1 


5 


8 


16 


19 


4 


— 


— 


35 (-5)3 


5 


9 


35 


35 


4 


0.01 s 


— 


72 (-2)'^ 


5 


10 


126 













5 


11 


126 


84 


10 


0.10 s 


0.03 s 


2« (-9)2 


5 


12 


246 


84 


92 


2.61 s 


0.26 s 


75 {-5y 


5 


13 


405 


90 


234 


14.22 s 


1.08 s 


103 (-3)1" 


5 


14 


550 


165 


44 


3.37 s 


0.44 s 


ll3 (-3)11 


5 


15 


1001 





996 




12.03 s 




5 


16 


1287 


78 


342 


2.43 m 


7.78 s 


313 (-3)13 


5 


17 


1820 





702 


5.48 m 


23.26 s 


lif (-7)1" 


5 


18 


2380 





3G4 


4.4() 111 


20.77 s 


14' (-4)'' 


5 


19 


3060 





192 


3.64 ni 


16.54 s 


154 (-4)15 


5 


20 


3876 





64 


2.13 m 


9.29 s 


791 191 (-l)i'' (-21)4 


5 


21 


4845 





64 


3.09 m 


13.53 s 


671 (-17)7 



Table 1: Data for branch-and-cut method with GLPK and CPLEX, for k e {3,4,5}. 
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7 

k 


n 


g[n,k) 


g{n,k) 


Nodes 


GLFK 




Stochastic 


Strong fcvxample 


6 


7 


1 


5 


2 




— 


— 


1^ (-6)1 


6 


8 


7 


14 


2 


— 


— 


— 


1^ (-7)1 


6 


9 


28 


28 


6 


0.01 s 


— 


0.03 s 


1** (-8)1 


6 


10 


70 


56 


32 


0.20 s 


0.05 s 


0.25 s 


8^ (-2)8 


6 


11 


126 


126 


10 


0.09 s 


0.03 s 


0.05 s 


9^ (-2f 


6 


12 


462 















6 


13 


462 


330 


24 


1.88 s 


0.16 s 


0.21 s 


2" (-11)^ 


6 


14 


924 


363 


294 


1.32 m 


6.17 s 


8.41 s 


212 (_i2)^ 


6 


15 


1705 


297 


12408 


2.94 h 


6.47 m 


4.76 m 


12^ (-3)1^ 


6 


16 


2431 


572 


1296 


34.19 m 


1.68 m 


54.20 s 


13^ (-3)1'^ 


6 


17 


3367 


1001 


266 


10.84 m 


44.14 s 


17.04 s 


143 (-3)14 


6 


18 


6188 





183960 




9.66 h 






6 


19 


8008 


560 


7262 




47.46 m 


13.74 m 


31^ (-16)^ 


6 


20 


11628 





27696 




4.58 h 


2.32 h 


^ 1 "7 — 7 ^ — \ 

3" (-17)^^ 


6 


21 


15504 





8932 




2.48 h 


1.11 h 


174 (_4)17 


6 


22 


20349 





4622 




2.06 h 


59.81 m 


18* (-4)i« 


6 


23 


26334 





2378 




1.66 h 


38.51 m 


194 (_4)19 


6 


24 


33649 





764 




49.35 m 




331 li« (-7)^ 


6 


25 


42504 





744 




1.15 h 


26.62 m 


IO4I 41^ (-21)** 



Table 2: Data for the branch-and-cut method using GLPK and CPLEX, for k = 6. 



k 


n 


g{n,k) 


g{n, k) 


Nodes 


Deterministic 


Stochastic 


Sharp Example 


7 


8 


1 


6 


2 






1' (-7)1 


7 


9 


8 


20 


2 






1« (-8)1 


7 


10 


35 


48 


8 






7^ (-3)'^ 


7 


11 


92 


118 


12 


0.02 s 


0.02 s 


3« (-8)^ 


7 


12 


246 


216 


100 


0.39 s 


0.92 s 


5'^ (-7)^^ 


7 


13 


462 


462 




0.36 s 


0.13 s 


11^ (-2)11 


7 


14 


1716 





26 








7 


15 


1716 


1287 


58 


6.02 s 


9.59 s 


21^ (-13)" 


7 


16 


3432 


1573 


1562 


4.97 m 


15.40 m 


214 (_i4)2 


7 


17 


6116 


1892 


26852 


2.61 h 


1.07 h 


125 (_5)12 


7 


18 


10296 


2080 


450772 


3.44 d 


1.02 h 


51^ (-13)^ 


7 


19 


14924 


3640 


28778 


1.25 d 


1.15 h 


163 (-3)16 


7 


20 


20944 


6188 


3615 


8.51 h 


28.80 m 


173 (-3)1'^ 


7 


21 


38760 













7 


22 


50388 


3876 


795236 


160.57 d 


3.08 h 


3i« (-19)^ 


7 


23 


74613 





13013* 




8.09 d* 




7 


24 


100947 





7870* 




4.61 d* 




7 


25 


134596 





6531* 




3.85 d* 




7 


26 


177100 





12718* 




9.25 d* 




7 


27 


230230 





30807* 




19.18 d* 




7 


28 


296010 





5564* 




2.51 d* 




7 


29 


376740 





6002* 




3.38 d* 





* These node counts and CPU times are averages of at least three runs using StochasticPropagation. 
Table 3: Completed computations for /c = 7 using CPLEX. 



k 


3 


3 


3 


4 


4 


4 


4 


5 


5 


5 


5 


5 


n 


11 


12 


13 


14 


15 


16 


17 


17 


18 


19 


20 


21 


g{n, k) 


45 


55 


66 


286 


364 


455 


560 


1820 


2380 


3060 


3876 


4845 


9s{n,k) 


46 


80 


84 


311 


375 


455 


750 


1946 


2562 


3165 


4876 


6097 



k 


6 


6 


6 


6 


6 


6 


n 


20 


21 


22 


23 


24 


25 


g{n,k) 


11628 


15504 


20349 


26334 


33649 


42054 


9s{n,k) 


12376 


17136 


21777 


27303 


39836 


50456 


Time 


5.71 d 


15.91 d 


2.26 d 


19.70 h 


67.60 d 


30.00 d 



Table 4: Comparisons of g{n^ k) and g's(n, k) when f{k) <n< f{k) + k. 



B A Computer-Generated Proof that /(3) < 11 

The fohowing proofs were created by executing BranchAndCut with propagation step Propa- 
GATePositive (Algorithm 3) and writing down the A;-sums which are determined to be strictly 
negative or nonnegative. 

Theorem 16. 3(8, 11) = (2°) = 45. 

Proof. The following sums generated by PropagateNegative (Algorithm 2) must be strictly 
negative or we have at least 45 of nonnegative sets: 

Xl + XQ+Xll X1+X8+X10 X2+X^ + Xll 

X2+XQ + XiQ X2+X7 + Xg X3+X4 + .Til 

X3+X5 + Xg X3+X7 + Xg X4+X6 + Xg 

The following sums generated by PropagatePositive (Algorithm 3) must be nonnegative or else 
the associated linear program (Definition 7) becomes infeasible: 

X4+X6 + X7 X4+X5 + X8 X3+X4 + Xio 

These positive sets now force at least 56 nonnegative 3-sums, and our target was 45 nonnegative 
3-sums. ■ 

Theorem 17. g{3, 13) = (2^) = 66. 

Proof. The following sums generated by PropagateNegative (Algorithm 2) must be strictly 
negative or we have at least 66 nonnegative sets: 

X1+X5 + Xi3 Xi+Xg + X12 X1+.X7 + Xii 

3^3+3^5 + Xi2 X3+X6 + Xio X4+X5 + Xii 

X4+X7 + Xg 

The associated linear program is infeasible. ■ 
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C A Computer-Generated Proof that /(4) < 14 

The following proofs were created by executing BranchAndCut with propagation step Propa- 
GATePositive (Algorithm 3) and writing down the A;-sums which are determined to be strictly 
negative or nonnegative. 

Theorem 18. g{i, 14) = (3^) = 286. 

Proof. The following sums generated by PropagateNegative (Algorithm 2) must be strictly 
negative or we have at least 286 nonnegative sets: 



Xi + Xj+Xis + xu 


x\ + a:8+xi2 + xi4 


X\ 


+ X9+X11 + Xi4 


X2 + X5+X10 + Xi4 


X2 + a^e+ajg + xi4 


X2 


+ X6+X10 + Xi3 


X2 + Xs+Xq + Xi3 


X3 + X4+X11 + a:i4 


X?, 


+ X5+X9 + Xi4 


X3 + X5+X10 + Xi3 


3^3 + a:6+-T8 + 3^14 


xz 


+ X6+X9 + Xi3 


X3 + Xq+Xio + X12 


x-i + xr+xs + xi3 


xz 


+ X7+X9 + X12 


X4 + X5+X8 + Xl4 


a;4 + xs+xg + X13 


X4 


+ X6+X8 + Xl3 


X4 + X6+Xg + X12 


X4 + X8+X10 + Xll 


Xf> 


+ X7+X8 + X12 


X5 + Xr+Xio + Xii 


xe + X8+X9 + Xll 







The following sums generated by PropagatePositive (Algorithm 3) must be nonnegative or 
else the associated linear program becomes infeasible: 



X3 + X4+X7 + Xg 
X\ + X7+X8 + Xg 
X\ + X5+X7 + X12 
Xi + X4+X8 + Xi3 
a^l + X2+X% + Xi4 



X3 + X4+X6 + Xio 
X\ + X6+X8 + Xll 
Xi + X4+X10 + Xll 
Xi + X4+X5 + Xi4 



X2 + X4+X8 + Xio 
X\ + X5+Xg + Xll 
Xi + X4+Xg + X12 
X\ + X2+X10 + Xi3 



These positive sets now force at least 199 nonnegative 4-sums. 

The following sums generated by PropagateNegative (Algorithm 2) must be strictly negative 



or we have at least 286 of nonnej 


native sets: 




Xi + X7+X12 + Xi4 


Xi + X8+X11 + Xi4 


X2 + X4+X11 + Xi4 


X2 + a;5+Xg + Xi4 


X2 + X^^Xxx + Xi3 


X2 + X^^X^ + Xi4 


X2 + X6+Xg + Xi3 


X2 + a;7+xio + X12 


X3 + X4+X10 + Xi4 


X3 + X4+X12 + Xl3 


X3 + X5+X8 + Xi4 


X3 + X5+X9 + Xl3 


X3 + X5+X10 + X12 


X3 + X6+X7 + Xl4 


X3 + X6+X8 + Xl3 


X-i + X6+X9 + X12 


X3 + X7+X8 + X12 


X3 + X7+X10 + Xll 


X2, + X8+X9 + Xll 


X4 + X5+X7 + Xi4 


X4 + X5+X8 + Xi3 


X4 + X5+X9 + X12 


X4 + X6+X7 + Xi3 


X4 + X6+X8 + X12 


X4 + X6+X9 + Xll 


X.f, + X7+X8 + Xll 
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The following sums generated by PropagatePositive (Algorithm 3) must be nonnegative or 
else the associated linear program becomes infeasible: 



X3 + X5+Xe + Xg 
X2 + X5+X8 + Xio 

xi + xr+xio + xii 

Xl + X5+X10 + Xl3 



X3 + X4+XS + Xio 
X2 + X4+Xg + Xio 
Xl + X7+X8 + Xi3 
Xl + X5+X8 + XU 



X2 + Xq+X-j + Xio 
Xl + X8+X9 + Xio 
Xl + X6+X9 + X12 
Xl + X4+X10 + Xl4 



These positive sets now force at least 265 nonnegative 4-sums. 

The following sums generated by PropagateNegative (Algorithm 2) must be strictly negative 
or we have at least 286 nonnegative sets: 



Xl + X5+X12 + Xi4 
Xl + X8+X10 + Xi4 
X2 + X3+X10 +X13 
X2 + X5+X6 + Xl4 
X3 + X4+X6 + Xl3 
X3 + X5+X8 + Xii 
X4 + X5+X6 + X12 
X5 + X7+X8 + Xg 



Xl + X6+X11 + Xi4 
Xl + X8+X11 + Xi3 
X2 + X4+X7 + Xi3 
X2 + X5+X8 + X12 
X3 + X4+X8 + X12 
X3 + X6+X7 + Xll 
X4 + X5+X7 + Xll 



Xl + X7+X12 + Xi3 

X2 + X3+X9 + Xi4 
X2 + X4+X9 + X12 
X2 + X6+X9 + Xll 
X3 + X5+X7 + X12 
X^ + X7+X9 + Xio 
X4 + X6+X8 + Xio 



The associated linear program is infeasible. ■ 
Theorem 19. ^(4, 15) = (3^) = 364. 

Proof. The following sums generated by PropagateNegative (Algorithm 2) must be strictly 
negative or we have at least 364 nonnegative sets: 



Xl + X8+X14 + Xi5 
X2 + X5+X11 + Xi5 
X2 + X7+X9 + Xi5 
X2 + X10+X11 +X13 
+ X5+X10 + Xi4 
X3 + X6+X10 + Xi3 
X4 + X5+X12 + Xi3 
X4 + X7+X10 + X12 



Xl + X9+X13 + Xi5 
X2 + X6+X10 + Xi5 
X2 + X7+X10 + Xi4 

X3 + X4+X13 + Xi5 
X3 + X6+X8 + Xi5 
X3 + X8+X9 + Xi3 
X4 + X6+X9 + Xl3 
X5 + X6+X8 + Xl4 



Xl + X11+X12 + Xi5 
X2 + X6+X11 + Xi4 
X2 + X8+X12 + Xi3 
X3 + X5+X9 + Xi5 
X3 + X6+X9 + Xi4 
X3 + X9+X11 + X12 
X4 + X7+X8 + Xi4 
X5 + X8+X9 + X12 



The following sums generated by PropagatePositive (Algorithm 3) must be nonnegative or 
else the associated linear program becomes infeasible: 



X3 + X4+X7 + Xio 
Xl + X7+X9 + X12 
Xl + X5+X8 + Xl4 
Xl + X3+X8 + Xi5 



X3 + X4+X6 + Xll 
Xl + X6+X8 + Xi3 
Xl + X4+X11 + Xi3 
Xl + X2+X11 + Xi5 



X2 + .X1+X7 + Xll 
Xl + X5+X11 + X12 
Xl + X4+X7 + Xi5 
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These positive sets now force at least 267 nonnegative 4-sums. 

The fohowing sums generated by PropagateNegative (Algorithm 2) must be strictly negative 
or we have 364 nonnegative sets: 



Xl 


+ X3+X13 + Xi5 


Xl 


+ X4+X9 + Xi5 


Xl 


+ X4+X11 + Xi4 


X\ 


+ X5+X8 + Xi5 


Xl 


+ X5+X9 + Xi4 


Xl 


+ X5+X10 + Xi3 


Xl 


+ Xq+Xs + Xi4 


Xl 


+ X6+X9 + Xi3 


Xl 


+ X7+X11 + X12 


Xl 


+ xs+xio + a::i2 


X2 


+ X3+X9 + Xi5 


X2 


+ X3+X10 + Xi4 


X2 


+ X4+X8 + Xi4 


X2 


+ X4+X9 + Xi3 


X2 


+ X5+X7 + Xi4 


X2 


+ X5+X8 + Xi3 


X2 


+ X5+X9 + X12 


X2 


+ xe+xg + X12 


X2 


+ Xg+Xio + Xii 


X2 


+ X7+X9 + Xll 


X3 


+ X4+X7 + Xi5 


X3 


+ X4+X10 + X12 


X3 


+ X5+X6 + Xi5 


X3 


+ X5+X7 + Xi3 


X3 


+ X5+X8 + Xl2 


X3 


+ X5+X10 + Xll 


xz 


+ X6+X8 + Xll 


X4 


+ X5+X9 + Xll 


Xa 


+ X6+X7 + X12 


X4 


+ X8+X9 + Xio 


X5 


+ X7+X9 + Xio 











The associated linear program is infeasible. ■ 
Theorem 20. g{4, 17) = (3^) = 560. 

Proof. The following sums generated by PropagateNegative (Algorithm 2) must be strictly 
negative or we have at least 560 nonnegative sets: 



Xl 


+ X6+X16 + Xi7 


Xl 


+ X7+X13 


+ 


Xi7 


Xl 


+ 


X8+X12 + Xi7 


Xl 


+ X8+X15 + X16 


Xl 


+ X9+X13 


+ 


Xi6 


X2 


+ 


a:5+a:i5 + a:i7 


X2 


+ X6+X12 + Xi7 


X2 


+ X6+X14 


+ 


Xl6 


X2 


+ 


X7+X11 + Xi7 


X2 


+ X7+X12 + X16 


X2 


+ X8+X10 


+ 


Xi7 


X2 


+ 


Xs+Xii + X16 


X2 


+ X8+X13 + Xl5 


X2 


+ X9+X12 


+ 


Xl5 


X3 


+ 


3:5+3^11 + a;i7 


X3 


+ X5+X13 + Xl6 


X3 


+ X6+X10 


+ 


xie 


X3 


+ 


X6+X12 + Xl5 


X3 


+ X7+X9 + Xi7 


X3 


+ X7+X11 


+ 


Xl5 


X3 


+ 


X7+X13 + Xi4 


X3 


+ X8+X10 + Xi5 


X3 


+ X8+X12 


+ 


Xi4 


X3 


+ 


X9+X11 + Xi4 


X4 


+ X5+X10 + Xi7 


X4 


+ X5+X11 


+ 


3^16 


X4 


+ 


X6+X9 + Xi7 


X4 


+ X6+X11 +X15 


X4 


+ X6+X13 


+ 


Xi4 


X4 


+ 


X7+X9 + X16 


X4 


+ X7+X10 + X15 


X4 


+ X7+X11 


+ 


Xi4 


X4 


+ 


X8+X9 + Xi5 


X4 


+ Xg+Xio + Xi4 


X4 


+ X9+X12 


+ 


a;i3 


X5 


+ 


X6+X9 + X16 


X5 


+ X6+X10 + Xi5 


X5 


+ X6+X12 


+ 


Xi4 


X5 


+ 


X7+X8 + Xi7 


X5 


+ X7+X9 + Xi5 


X5 


+ X7+X10 


+ 


Xi4 


X5 


+ 


Xs+Xii + Xi3 


Xq 


+ X7+X12 + Xi3 


X7 


+ X9+X10 


+ 


a;i3 









The following sums generated by PropagatePositive (Algorithm 3) must be nonnegative or 
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else the associated linear becomes infeasible: 



X4 + XQ+Xg + Xii 
X4 + X5+Xg + Xi4 
X3 + X6+X9 + Xi3 
X3 + X5+X10 + Xu 
X3 + X4+X9 + Xi5 

These positive sets now force at 
4-sums. 



X4 + Xq+Xj + ,Ti3 
X4 + X5+X6 + Xi5 

a:;3 + xe+X7 + 2:14 

X3 + X5+X7 + Xi5 

X3 + X4+xe + xie 
560 nonnegative 4-sums, 



X4 + X5+a;io + xii 

X3 + + Xll 

X3 + Xs+Xii + Xi2 
X3 + X4+X11 + Xi3 
X2 + iCs+Xio + X12 

our target was 560 nonnegative 
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